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Целью работы является изучение новых математических мо- 
делей, включающих в себя генератор колебаний и полуогра- 
ниченное неоднородное по глубине основание, обладающее 
пористостью, флюидонасыщенностью, вязкоупругостью. В 
качестве основания рассмотрены пористоупругий слой, 
насыщенный смесью жидкости и газа, гетерогенный слой с 
вязкоупругим покрытием, гетерогенный слой с заглублен- 
ным жидким слоем. Основание пакета слоев — жесткое. 
Действие поверхностного осциллятора представляется в 
виде ряда Фурье и решается задача с установившимся ре- 
жимом колебаний. Применение интегрального преобразова- 
ния Фурье к уравнениям, описывающим сплошные среды, 
при удовлетворении граничным условиям позволило по- 
строить интегральные формулы, описывающие напряженно- 
деформированное состояние в пакете слоев. Предложен чис- 
ленный алгоритм для изучения зависимости распространения 
поверхностных волн от механических и геометрических ха- 
рактеристик задачи. Описанные модели широко применяются 
в геофизике, сейсморазведке, строительстве, проектирова- 
нии железнодорожных магистралей, конструировании но- 
вых материалов. 


Ключевые слова: гетерогенная слоистая среда, волновое 
поле, распространение колебаний, заглубленный жидкий слой. 
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Введение. Моделирование сложных технических динамических систем, содержащих сосредоточенные и 
непрерывно распределенные параметры, представляет научный интерес с точки зрения использования в геологии, 
сейсморазведке, строительстве, проектировании железнодорожных и автомагистралей. Применение таких моделей 
позволяет совершенствовать биотехнологии, конструировать новые материалы с заданными свойствами. Следует 
отметить, что требования к проектированию, эксплуатации сложных технических объектов и технологических про- 
цессов растут. В связи с этим возникает необходимость изучения новых математических моделей динамических про- 
цессов, описываемых системой обыкновенных дифференциальных уравнений, моделирующих воздействие генера- 
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тора колебаний, и дифференциальных уравнений в частных производных, описывающих поведение полуограни- 
ченного основания. В качестве основания в настоящей работе рассматриваются случаи гетерогенного слоя (задача 
А), гетерогенного слоя с вязкоупругим покрытием (задача В), гетерогенного слоя с заглубленным жидким слоем 
(задача С). 

Основная часть. 

1. Рассматривается плоская задача о колебаниях гетерогенного слоя (задача А), гетерогенного слоя с вязко- 
упругим покрытием (задача В), гетерогенного слоя с заглубленным жидким слоем (задача С). Полагаем, что дей- 
ствие генератора колебаний представлено в виде ряда Фурье. Поэтому к лицевой поверхности составного слоя при- 


—Г®Е 


ложена нагрузка Р(хуе ‚—а <х<а, у=й , осциллирующая с частотой о. 








— | 0, х >а, и 
а(х, у.) |= 
№ Ве“. х<а 
где 4(х, у,1) ={0.у,б„,} — вектор напряжений в верхнем слое. 


В верхнем гетерогенном слое нагрузка приложена по типу «непроницаемый поршень», нижний слой лежит на 
недеформируемом основании. 
Вязкоупругое покрытие слоя описывается уравнениями Ламе [1], коэффициенты которого являются величи- 
нами с малой комплексной частью в соответствии с гипотезой внутреннего трения Сорокина: 
д^и, 
д’ 


Для учета пористости среды, ее Газонасыщенности и водонасыщенности чаще всего применяется модель Био 





ии, +(^+и),;=р ; 1=Ь2. (2) 


— Френкеля [2, 3], уравнение перемещения для которой имеет вид: 

















и; Аб ди ди |, 
Р11 РУ Р12 РУ 5, 5 0. 
д?и. ду, ди. ОУ 
+ —"-_Ь -— =с7, ;1=12, 3 
Р12 5 Р22 52 и 2 (3) 


где и; (х, у, <,1),у, (х, у, 5,1) 1=1,2 — компоненты векторов перемещений упругого скелета и флюида, р!› <0 — коэф- 


фициент динамической связи упругого скелета и жидкости. 
Коэффициенты динамической плотности выражаются формулами: 


ри =И-т)р, — р12, 
Р>2 = 1—т)р ; —Р12, 


2 
т 
где р;, р; — плотности сред упругого скелета и жидкости, т — пористость среды, = —. | — коэффициент 


О 


вязкости жидкости, К, — коэффициент проницаемости. 


Смешанная краевая задача для пористой среды другой модели рассматривалась в [4]. 


Связь между тензором полных напряжений Г, и деформациями упругой е, и жидкой =, фаз выражается в 


виде: 
о" = Аеб, +2№, +08, с/=О+КЕ, 
9 = ли, В = ай», (4) 
Гу =о, +60 7, 5 1=1,2, 
где д, — символ Кронекера, д, =1,6,, =01=; в, — тензор напряжений, действующих в упругом скелете; с/ — 
напряжения, действующие на жидкость в порах; А, М, О, К — механические характеристики гетерогенной среды, за- 


висящие от скоростей распространения волн в упругом скелете и в жидкости. 


Вектор перемещений ии, (х, у,1)} и давление в идеальной жидкости заглубленного слоя выражаются через 


волновой потенциал ф (х, у,ЁК): 


ОФ (х,у,1). 


Ро(х, у,Ю) = -Ро т их, у,1) = атаа ф (х,у,1). (5) 
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Так как режим колебаний предполагается установившимся, отделим временной множитель 


> —_* — а 
и(х, у) =и (х, у)е '°’. В дальнейшем опускаем знак * и рассматриваем амплитудные значения соответствующих 


величин. 
В задаче А на нижней границе слоя равны нулю нормальное перемещение и касательные напряжения, а также 
выполнено условие непроницаемости нижней границы: 


и (х, У) | = У, (х, у) | = 0; 


у=-й> 





из (У) [и => (%У) | — =0:[х| <а, (6) 


9906) | =0; #1. 
В задаче В на границах раздела слоев задаются условия сцепления, для гетерогенного и вязкоупругого слоев с 


номерами К, К +1 соответственно граничные условия имеют вид: 


—К+1 
+ 


и" (ху) | =и (№) |; 


К К 
из (х,У) в =2 (№5) в; 


К К+1 . 
©1> (х, У) реа = 612 (ху) | и 
у. т. (7) 
(652 (х,у)+01 (х,У)) | и, = 


К-+1 
=0> (х, у) |. 


В задаче С на границе заглубленного жидкого слоя и пористоупругой среды предполагается свободная филь- 
трация жидкости через границу. При этом требуется: 


(1-ти (х, У) |6 + ту (х, У) |5 = и» (х, У) ре: › 
(055 (х,у)+055 (х,5)) |0 =(т-Юрь | о, (8) 


‚К р ‚К 
025 (%,У)) |о=-тро [40$ 915 (4%) [и =0. 


2. Волновые поля поставленных задач находятся применением преобразования Фурье: 
— оо . 
((а, у) = | и(х, уде" "ах. 
со 

Представления вектора перемещений в виде вихревой и потенциальной части и использование свойств преоб- 
разования Фурье и дифференциальных операторов позволяют свести решение систем уравнений (1)-—(8) к волновым 
уравнениям. Для описания упругой среды достаточно двух потенциалов, для пористоупругой необходимы 3 потенци- 
ала среды. Подробнее эти преобразования описаны в работах [5], [6]. 

Определяя произвольные постоянные из граничных условий, приходим к описанию перемещений в составной 
полосе. Однако полученные соотношения должны быть преобразованы к виду, подходящему для дальнейшего реше- 
ния задач о составной полосе, не иметь неопределенностей и обеспечивать устойчивый счет для достаточно больших 
значений аргумента. Описанный алгоритм громоздких преобразований реализован в пакете Мар/е с помощью опера- 
ций матричной алгебры Пипа: 


| === _ 
и(х, у) = РЕ ГС (а, у) О(ае “аа; (9) 
ы 


О(а, у) = т уе“ ах . 


Элементы подынтегральных матриц-функций С (ау) в представлении (9) являются осциллирующими, убы- 


вают на бесконечности степенным образом, в комплексной плоскости — мероморфные, имеют вблизи вещественной 
оси счетное количество полюсов с малой мнимой частью. 


Контур интегрирования У, выбирается в соответствии с принципом излучения [5], обходя регулярные по- 


ложительные особенности подынтегральной функции в нижней комплексной полуплоскости. 
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Ау (а) 


Ввиду громоздкости элементов матрицы С (а, у) приведем значение для элемента С. (а) =. ( ‚ необхо- 
(И 





димое для вычисления вертикальных перемещений составной полосы в случае нормально приложенной нагрузки. 
Задача А (гетерогенный слой): 


п? Г 


5% = о -к”;: к ее (10) 
У, 
А(0) /2 = (713(8152611512 — 8251612511) + 2515253433 (ти — Оси — 
—(т, — е58н)@ св + (155 (т› -П)- 851 (т, — 1))54С1 112513; 


До (9) = (т — т, кз 51559 Р(0); 


8; = 24335; — КТ. Тно = 41 +4 — 243 + (42 +4»)т; #=1,2; 
51; =1—ехр(—25,И); с; =1+ехр(-—25,П); 
Чи =(А+2№)/Н,4›О!Н, 4 =К!Н,4в =2М№!Н, 
Н =А+2М№М+О+КБ, 54 =5 ва, 13 = 1+2 / 722. 
Здесь У, — скорости распространения волн трех типов в гетерогенной среде. 


Следует отметить, что при Й —> с, 51; Ес; =| ‚ формулы (10) для слоя вырождаются в формулы для гетероген- 
ного полупространства с непроницаемой границей. 
Задача В (двухслойная среда): 
А(а) = 463 +4163 +45, Ао(0) = (4563 + 4.53 +4, )Р(а), 
Ар = (61, +и° 5651 )6 15°, 
4 = [(ез6› +0166 | + 6>64Б>> )01 + (636 +0>636>› +61 64В 1 6 |616>, 
45 =е55 +[-61 65 (5 +5) +2а-ЬЬ, |616 +е.5', 
е. = а”, о 05 ы=Ь2, е. = а” (В. -Ь. |), 
Ч. =—01е4 (Ь:6> +Ь.01 6162, а =45, (ще 65 +Б:е>9} ), 
4; =Ь,>(е/55 +е.6\) +6164, (11) 
ед = 65 (65° +5” )610> +4564 1 (©1 )* +29” Ь,Б>> + 2о1ез (655 +64), В = (2453514 + 13 (Чл 442) 9 +(1-Ь)53, 
В = 2((5, 51 Уз +1 о +26 — Юз, В =-ПУз - 26433 +1 -6;6>> =В, 
03 =(Д-Р {ег + 53) — 2171з —4В, 
1 = 243347 (52—51 )— 94152 + 44291 › 
Ь = 243307 (52то — $1 ) — 44152т> + 442 Ту, 
[5 = 515253 (т; — т. ), 4 = (ть —1)5> —(т -15, 
Ь = 2(ус <> -61665с1), 5 =2(усьср 616541 со), 
6; =05 ас, 4 =6, 
91 =7° сс: —и*о1 02 си, бо = у сос фи’ 61 6с с, 
2 


+ 25 (2 2 — | Р.ь В О ее 
С в == ° у=и +о> Ох = Ц/( — 0х; 0. пу. Е а 
| 


Здесь У, — скорости продольных и поперечных волн. 


Задача С (гетерогенный слой с заглубленным слоем жидкости): 
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Д(а) = ЯР: —- 272 +В; 
Ар (0) = (52559 +5272 +48 + (из + вс 1+ 
0400) = (51 7251 т 52 74255 т 14453 ТАЙ т 75 )С] 1 
2 
(д — Из )со1 — 7 / Уб®осз1 )Р(а); 
= : = 2. : —. в , 
т = 757 — 6 +1; 752 = ЩИ +0 (ВПО +п1); 3 = ПВ +9 (5по — И); 
— 2 2 ТИВ р В 2 . 
Уз+: = 965; — К; (44 +45т‚)1=12; аб; = -К; (412 + 42от;}; 
2 
гороткзсШТуой (1— т-+ тт, ) 5 22 
= № ла -Ю,.Ко =©/\; 
274Н 
уз = (975 — 484) / (4648 )> о = (4 туб ут) / (тав), 5; = 5; / 5;; 
44 = 411 +Чр — 96;45 = 412 + 422 .ти =1-т;; (12) 
ба , №8 м — , 
п =О (51 51 + 55 5240 ) + 76353» № = УС + 4975 ро + 464813; 
—_ р О = ПЕ 2> / __ 2 > 
в, = "11151 51 + 51525249 + Тз1Уз53 > 1 = 9 (У751 52 /Уб -Сз) +70 53; 
Е , 
5 = (7577С12 + 46\оаз — 9 53.4683) /Уб; 
_ - - 2 
5 = (то у2755 + ть 983 — ТС) / Ув; 
7 = Т-1С12 Тс, В = 7451 — 7552, = ТТ 451612 — Ту 552611, 


По = 2 +СиьП1 = (И — ть соси, По = 491612 — 25261, Пз = И — 7; 


— 2  _ | 
Роза: = ПЗИП — 53, = 12; 
__ Т-+1 . 
по = (-0” 4ов; +77 /Убв;; 


По: = (Оу, + То / Уб:. 


Здесь У, — скорость распространения волн в жидкости. 


Подынтегральные функции (10)—(12) представлены как отношения целых функций в виде, обеспечивающем 
устойчивый счет интегралов на полубесконечном промежутке. Нули дисперсионных функций — знаменателей 
подынтегральных функций определяют скорость волн рэлеевского типа, распространяющихся по лицевой поверхно- 
сти среды. 

3. В ближней зоне интегралы вычисляются интегрированием по комплексному контуру. В области, дальней 
от приложения нагрузки, соотношение (9) можно вычислить по вычетам [7]. При этом значения полюсов находится 


численно из решения дисперсионного уравнения Л(9)=0. Первый полюс порождает распространяющуюся по по- 


верхности волну, которая переносит наибольшую долю энергии. Величина полюса определяет скорость распростра- 
нения волны. 
Численный анализ проводился при следующих механических характеристиках, определяющих вязкоупругую 


и гетерогенную среду [8—10], что соответствует водонасыщенному песчанику: 


р=1700 кг / м”, И, =120м/с, И, =208 м/с; 
й, = 0,4 м, й, =0,1 м, а=1 м; 


р, =2100 кг/м’, И =160 м/с, Г, =277 м/с; 
О=0,5585.10’ па, В =0,1746. 107 па; 
А=0,4272.10° па, № =0,4311.107 па, р; =1000 кг / м. 


Для всех типов оснований численный анализ дисперсионных соотношений выявил определенные закономер- 
ности. При увеличении жесткости скелета обводненного основания скорость распространения поверхностных волн 
увеличивается. Увеличение пористости гетерогенной среды, увеличение плотности флюида, насыщающего поры, 
напротив, приводят к уменьшению скорости распространения волн по поверхности среды. Эти закономерности иллю- 


стрирует рис. 1. 
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Рис. 1. Изменение полюса двухслойного основания при изменении жесткости подстилающего гетерогенного 
полупространства У5 и плотности жидкости в порах го } 


Наличие жидкого слоя значительно увеличивает динамичность среды, особенно при низких частотах. На рис. 
2 представлены графики, иллюстрирующие распространение волн по поверхности гетерогенного слоя с заглубленным 
жидким слоем — верхние кривые, и гетерогенного слоя на жестком основании — нижние кривые. По вертикальной 
оси отложены безразмерные перемещения, по горизонтальной оси — безразмерное расстояние от области приложения 


нагрузки. 
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Рис. 2. Зависимость действительной (сплошная линия) и мнимой (штриховая линия) части вертикальных перемещений 
поверхности среды и от расстояния х до области приложения нагрузки 


Описанный метод решения задач применим и к соответствующим задачам в пространственной постановке. 

Выводы. При изменении геометрических и механических характеристик любого из составляющих слоев ос- 
нования, обладающего пористостью, флюидонасыщенностью, вязкоупругими свойствами, изменяется волновое поле 
на поверхности пакета слоев. Применение простейших моделей оснований типа полупространства Винклера приводит 
к искажению качественной картины динамического процесса и не может использоваться при моделировании проблем 
динамики. 

Авторы выражают благодарность профессору ЮФУ М. А. Сумбатяну за внимание к работе. 
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